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ABSTRACT 

In a previous paper we presented measurements of the properties of jets and cores in 
a large sample of FRII radio galaxies with z < 0.3. Here we test, by means of Monte 
Carlo simulations, the consistency of those data with models in which the prominences 
of cores and jets are determined by relativistic beaming. We conclude that relativistic 
beaming is needed to explain the relationships between core and jet prominences, and 
that speeds between 0.5 and 0.7c on kpc scales provide the best fits to the data. 
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1 INTRODUCTION 

There is strong and widely accepted evidence for relativistic 
bulk speeds on parsec scales in radio galaxies and quasars. 
In the most extreme objects the observed superluminal 
motions, the rapid variability and consequent high bright- 
ness temperatures, and the absence of very strong inverse- 
Compton emission in the X-ray seem to require relativistic 
bulk speeds with 7^2; it seems reasonable to assume that 
such high bulk Lorentz factors are present close to the active 
nucleus in all extragalactic radio sources. Evidence that rel- 
ativistic speeds persist on kiloparsec scales is less conclusive, 
but in FRII objects a one-sided jet on parsec scales usually 
connects to a one-sided jet on kiloparsec scales, and the re- 
sult of Laing (1988) and Garrington et al. (1988) that the 
brighter jet in an FRII radio source tends to lie in the less 
depolarized lobe is very hard to explain unless relativistic 
beaming is still important on these scales. 

Relativistic beaming is an essential ingredient in unified 
models for powerful radio galaxies and quasars (e.g. Barthel 
1989); on the assumption that quasars are radio galaxies 
oriented at less than some critical angle to the line of sight, 
relativistic bulk speeds in the core and jet imply that the 
cores and jets of quasars should be more prominent rela- 
tive to the extended emission and that their jets should be 
more one-sided. Until recently, the lack of detected jets in 
FRII radio galaxies, while giving qualitative support to this 
model, meant that it could not be tested quantitatively. 

In a previous paper (Hardcastle et al. 1998; paper I) we 
discussed the properties of radio jets and cores in a large 
sample of SCR FRII radio galaxies with z < 0.3 which had 
mostly been imaged with the NRAO VLA at high resolu- 
tion and sensitivity. Since there are no quasars in the parent 
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samples at this low redshift, there is necessarily a popula- 
tion of galaxies in our sample aligned at small angles to the 
line of sight. We suggested that our results were consistent 
with a modified version of unified models in which broad- 
line radio galaxies (BLRG) were the aligned objects and the 
low-power counterparts of quasars; the jets and cores of the 
BLRG were systematically more prominent than those of 
narrow- line radio galaxies (NLRG), and the relative num- 
bers of NLRG and BLRG were consistent with their being 
divided by a critical angle to the line of sight of ~ 45°, 
similar to that found by Barthel (1989) from the relative 
numbers of higher-redshift SCR radio galaxies and quasars. 
The radio properties of the BLRG in our sample were sim- 
ilar to those of the FRII quasars studied by Bridle et al. 
(1994). Assuming that relativistic beaming is the dominant 
factor determining the prominences of the cores, the rela- 
tionships between the properties of cores and jets implied 
that the jets were beamed, with characteristic speeds in the 
jets lower than those in the cores, consistent with the re- 
sults of Bridle et al. (1994). In this paper we attempt to 
make this conclusion quantitative by fitting empirical mod- 
els of relativistic beaming effects to the data derived from 
the sample. 

The observed flux density Sohs of a jet with speed v = 
/3c and Lorentz factor 7 = (1-/3^ )~2^ emitting isotropically 
in its rest frame, is related to the flux density that would be 
observed if it were at rest (^rcst) by the relation 

-/3cose)]-('"+"> 
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(Ryle and Longair 1967) where 9 is the angle made by the 
velocity vector with the line of sight, a is the spectral index 
(defined throughout in the sense S oc and m is a 

constant reflecting the geometry of the beamed component 
(e.g. Scheuer & Readhead 1979); we use m — 2, the value 
appropriate for a continuous jet. Because in practice the 
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material in beams is expected to have a range of speeds, care 
is necessary in interpreting values of 7 and 13 deduced from 
this expression as straightforward physical quantities; they 
nevertheless provide a means of characterising the physical 
parameters in the beamQ 

The approach to modelling relativistic beaming that 
we adopt has been widely used in testing unified models of 
different classes of radio sources (e.g. Scheuer & Readhead 
1979; Urry & Schafer 1984; Urry, Padovani & Stickel 1991; 
Morganti et al. 1995). The prominence of a radio component 
is defined as the ratio of the flux density of the component to 
the flux density of the extended emission (which we define as 
the total fiux density of the source minus any detected emis- 
sion from the core and jets). It is assumed that the extended 
emission is unaffected by beaming; the prominence p' of a 
beamed component is then taken to be given by the prod- 
uct of some intrinsic prominence p and the Doppler beaming 
factor given by equation |l| 

p'=p[7(l-/3cose)]^("'+"' (2) 

To see whether observed data are consistent with a particu- 
lar set of assumptions about the parameters of such a model, 
we generate large numbers of simulated sources by Monte 
Carlo methods (assuming an isotropic distribution of angles 
to the line of sight) and compare the resulting distributions 
with the observations. By varying the unknown parameters 
(p and P in this simple model) and examining the goodness 
of the fit, it is possible to find the values of the parameters 
that produce simulated datasets that, on average, best re- 
produce the data (i.e., have the highest probability, on some 
statistical test, of being drawn from the same parent popula- 
tion). Depending on the fitting statistic used and the space 
over which it is maximised, this approach may or may not 
be formally identical to a maximum-likelihood analysis; but 
in any event we expect it to provide insights into the regions 
of parameter space suggested by the data. 

In the rest of the paper we apply this simple technique 
to the data of paper I. 



2 THE DATA 

The sample of paper I consisted of 50 radio galaxies from 
the samples of Laing, Riley & Longair (1983) and the SCR 
catalogue (Spinrad et al. 1991). All had redshifts less than 
0.3 and radio power at 178 MHz greater than 1.5 x lO^'^ W 
Hz~^ sr~^, above the canonical FRI-FRII break (Fanaroff 
& Riley 1974) although a few had structures which resem- 
bled those found in FRIs (e.g. apparently dissipative jets). 
The reader is referred to paper I for details of the selection 
criteria. In what follows we assume, as argued in that pa- 
per, that the selection was not strongly biased in ways that 
might affect unified models; in particular, that objects at all 
angles to the line of sight are represented in the sample. 

High-resolution radio maps were available to us for 44 
of the 50 objects, and from them we were able to measure 

t As in paper I, we try to follow Hughes & Miller (1991) in dis- 
tinguishing between the 'beam', the stream of particles supplying 
the hot spots and lobes with energy, and the 'jet', its observable 
manifestation. 



the flux densities of cores and jets and counterjets, or to set 
upper limits on these quantities where the components were 
not detected. These quantities, are reproduced from paper 
I in Table |l| As argued by Bridle et al. (1994), the straight 
component of the jet is the only part in which relativistic 
beaming effects can usefully be studied, and throughout this 
paper 'jet flux density' and 'jet prominence' will be used to 
refer to the flux density and prominence of the straight com- 
ponent of the (brighter) jet only, as defined in paper I. We 
do not expect to introduce any bias into the data by con- 
sidering the straight jet only. The prominences of the core 
and straight jet are also tabulated. If more than one jet is 
detected in a given source, the prominence of the brighter 
jet is used; if only one jet is detected, we tabulate its promi- 
nence; if no jets are detected, the larger of the two upper 
limits on jet fiux density is used to derive an upper limit on 
jet prominence. The errors quoted for the fluxes of detected 
jets are intended to reflect the uncertainty in the separation 
of the jet from its background, as discussed in paper I, and 
are not formal errors. Jets or possible jets were detected in 
33 of these 44 objects; radio cores were detected in all but 
2. 

We also tabulate the optical emission-line classes of the 
sources. We argued in paper I that there is a real difference, 
in this sample, between the properties of those sources with 
low-excitation emission lines (low-excitation radio galaxies, 
LERG) and those with stronger emission lines (the narrow- 
line radio galaxies, NLRG, and broad-line radio galaxies, 
BLRG). The prominences of cores and Jets for the different 
classes of object are plotted in Figures ul and H Fig. I shows 
that the most prominent jets are actually those of a sub- 
set of the LERG; in paper I we suggested that these were 
sources with particularly dissipative beams. For this reason 
we restrict ourselves in what follows to the sub-sample of 31 
high-excitation objects. 



3 MODEL FITS 

Both the distributions of core and jet prominence (plotted 
in Figs |l] and ^ and the relationship between these quanti- 
ties (Fig. ^ contain information about the best-flt beaming 
parameters. It is clear that the distributions could individ- 
ually be consistent with a relativistic beaming model while 
exhibiting some relationship (for example, a core-jet promi- 
nence anticorrelation) that effectively ruled such a model 
out. In paper I we showed that both the distributions of 
these quantities and the positive correlation between them 
were qualitatively consistent with a model in which both 
cores and jets are relativistically beamed, and with unified 
models in which the BLRG were the objects more closely 
aligned with the line of sight. Traditionally the procedure 
outlined in Section |l] has been applied to a single distribu- 
tion, e.g. to core prominences (Morganti et al. 1995) or jet 
sidednesses (Wardle & Aaron 1997). We begin our quantita- 
tive approach to the problem in that way; later we discuss 
the implications of the core-jet relationship. 

A serious problem in the data analysis arises because of 
the number of upper limits in the distributions, particularly 
in the jet prominences. Whether a jet is apparent in a par- 
ticular map or not depends not just on the sensitivity of the 
observations, but also on the presence or absence of confus- 
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ing features in the lobe; it is therefore difficult to reproduce 
the non-detections in the data in numerical simulation. We 
believe that the upper limits we determine are probably not 
far above the true values, but in what follows we shall try 
to assess the effects of the upper limits on the results of the 
fits. 

In order to apply the model discussed in section |l| to 
a large number of sources, it is necessary to make some as- 
sumptions about the intrinsic prominence p of the compo- 
nent of interest. The simplest assumption, adopted by Urry 
& Schafer (1984) and subsequently followed by other au- 
thors, is that the intrinsic prominence of jets or cores is 
fixed (equivalently, that the intrinsic luminosity of jets or 
cores is a fixed fraction of the total luminosity of the radio 
source) so that all sources have the same intrinsic promi- 
nence for a given component. This is a reasonable assump- 
tion so long as relativistic beaming is expected to be the 
dominant effect in determining the prominence of a compo- 
nent. In the case of radio galaxies, it seems likely that the 
environment of, and 'weather' in, a particular radio source 
can strongly affect the efficiency of jets and cores, and it is 
certainly true that the environment will affect the extended 
luminosity of a source. The normalisation by extended flux 
density removes much of the dependence of the prominences 
of jet and core on the underlying power of the radio source 
('beam luminosity') but cannot eliminate it entirely. All of 
these factors affect component prominences, in general by 
introducing additional scatter into the distribution of ob- 
served prominences. In our analysis we begin by making the 
simplest assumptions possible, and go on to discuss the ef- 
fects of relaxing them. 

The statistical test we use for comparing real and simu- 
lated data throughout the paper is the Kolmogorov-Smirnov 
(K-S) two-sample test; the statistic provided by this test can 
be converted to a probability that the two samples compared 
are drawn from the same parent population [see e.g. Press 
et al. (1992) for a discussion of the properties and imple- 
mentation of this test]. We prefer this test to the test 
because of the loss of power involved in binning our small 
sample; however, wherever we compare two one-dimensional 
distributions in the following sections we find that the re- 
sults when a test is used are consistent with those of the 
K-S test. 

As discussed in section |2|, we believe that the data are 
not biased with respect to their orientation to the line of 
sight; this allows us to neglect the details of unified models 
and to treat the 31 NLRG and BLRG as a homogeneous 
sample. The clear differences, discussed in paper I, between 
the core and jet prominences of the BLRG and NLRG moti- 
vate the belief that relativistic velocities are involved on both 
scales; but it would also in principle be possible to use the 
presence or absence of broad lines as an additional source of 
information on the inclination of sources in the sample when 
carrying out model fits. We have chosen not to do this for 
two reasons: firstly because the presence or absence of broad 
lines in these objects is strongly dependent on the quality 
of available optical observations (see e.g. Laing et al. 1994) 
and because this raises the question of the status of objects 
in the sample, like 3C 234, with polarized broad lines; sec- 
ondly because to generate comparison samples requires a 
knowledge of the critical angle separating the two popula- 
tions, which, even if such a simple unified model adequately 



represents the true situation, is not well constrained, given 
the small sample size, by the existing data on numbers of 
sources, supposing the numbers themselves to be accurate. 
We therefore restrict ourselves, in this section, to noting the 
degree of agreement between the data and the predictions 
of a simple unified model (with critical angle 45°) based on 
our best-fit models. 

Throughout the rest of the paper a subscript c denotes 
a quantity related to the core of a source and a subscript j 
denotes one related to the jet. 

3.1 Cores 

Only one of the core prominences for the high-excitation 
objects is an upper limit (that for 3C 153) and therefore we 
ignore the problem of upper limits in this part of the analy- 
sis, treating the upper limit on 3C 153's core prominence as 
a measurement. We assume that the spectral index of cores 
is zero. Because cores are the unresolved bases of jets, which 
we presume to be intrinsically two-sided, the correct model 
to use for our sample, where there are angles to the line of 
sight 6 close to 90°, is 

Pc = { [7(1 - Pc COS e)]-('"+"' + [7(1 + I3c cos e)]-('"+-) } (3) 

— this is identical to equation ^ at small 6 and large 
/3c, where the second term is negligible in comparison to the 
first. 

Fig. ^ shows the distribution of core prominences in the 
NLRG and BLRG to be very broad. It is therefore not sur- 
prising that when models of the form described in section 
|l| using the prominence described by equation ^, are fitted 
to the data only those with large values of /3c can reproduce 
the data well (Fig. ^). The data do not put an upper con- 
straint on /3c, but they are best reproduced with /3c > 0.95 
(corresponding to 7c ^ 3). The lack of an upper constraint 
is due to the relatively few highly beamed objects in a sam- 
ple without orientation bias. VLBI information constraining 
the speeds in the cores is only available for two of our ob- 
jects: superluminal motion in the core of 3C 111 constrains 
7 6.9 (Preuss et al. 1990) and in 3C 390.3 7 3.6 (Alef 
et al. 1996). Our best-fit Lorentz factors are consistent with 
these observations. For /3c = 0.98, under our assumed unified 
model with a critical angle of 45° , we would expect the me- 
dian BLRG core prominence in a sample of 31 objects to be 
about 18 times the median NLRG core prominence, which is 
larger than the observed ratio of about 8; however, the ratio 
observed in our data is not inconsistent with this value of 
the critical angle at the 95 per cent confidence level for any 
value of /3c allowed by the distribution of core prominences, 
though it is definitely inconsistent with a critical angle ^ 40° 
unless /3c < 0.98. 

In paper I we discussed a possible positive correla- 
tion, only apparent in the narrow-line objects, between 
core prominence and apparent source linear size. Regression 
[performed using the Buckley- James method, implemented 
in the survival analysis software package ASURV Rev. 1.1; 
LaValley, Isobe & Feigelson (1992)] suggests a non-linear re- 
lationship between the two quantities, p'^ oc (L sin ^)^ ''*'''*. 
This might be attributed to a proportionality between un- 
projected source length L and intrinsic core prominence pc; 
such an effect would tend to be seen preferentially in the 
narrow-line objects, since it would be diluted by projection 
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and beaming effects in objects at small angles to the line of 
sight. We may consider two hypotheses: 

(i) Intrinsic core prominence is independent of source 
length; prominences only depend on /3c and 9 as in equa- 
tion ^ and the observed correlation is coincidental. 

(ii) Intrinsic core prominence is dependent on source 
length; prominences depend on (5c, 9 and L, the unprojected 
source length. 

In the latter case the core prominence of a source is given 
by 

= fc,,iL" {[7(1 - 13c cos + [7(1 + Pc cos } 

where a governs the relationship between L and p'c- The 
effects of projection and beaming will tend to flatten the 
power-law slope of the observed relationship compared to 
the true one. The magnitude of this effect, observed in sim- 
ulated datasets, depends on /3c and the unification angle, but 
for plausible values of both we expect the observed slope to 
be flatter than the true slope by about 0.2, so that the true 
slope a is likely to be about 1.8 ± 0.4. For simplicity, there- 
fore, we choose a = 2, a value consistent with the data. fcc,i 
is the constant of proportionality in this relationship and 
clearly (for this choice of a) has units of kpc~^. 

To use Monte Carlo methods to test these models 
against the data, we first need to model the distribution of 
unprojected source lengths in the parent population. Lit- 
tle is known about this distribution. The simple empiri- 
cal form suggested by Kapahi (1976) is not a particularly 
good fit to our data, as it overpredicts the number of small 
sources. Simulation shows that the data are moderately well 
described by a simple normal distribution of unprojected 
length (mean 560 kpc, standard deviation 230 kpc) with neg- 
ative lengths being discarded. We use simulated unprojected 
lengths drawn from this distribution to test the hypotheses. 

We first repeated the core prominence simulations to 
find the best-fit values of A;c,i and /3c for model (ii). The re- 
sults of the simulations for /3c are not strongly affected by 
the inclusion of a length dependence, presumably because 
the distribution of core prominences in the data is much 
broader than that of lengths. Lower values of /3c become 
slightly more probable, as we would expect (Fig. . There 
is again no constraint on the maximum value of the beaming 
speed /3c, so we choose a representative value of /3c = 0.98, 
corresponding to fcc,i ~ 2.5 x 10~^ kpc~^. Using these values, 
we were able to generate pairs of (projected length, apparent 
core prominence) values to match those in the real dataset. 
We can then answer the question: how likely is it under the 
two hypotheses that 22 NLRG cores would produce the sig- 
nificant positive correlation seen on a Spearman Rank test 
in the real data? Assuming that a critical angle of 45° to the 
line of sight divides broad- and narrow-line objects, and set- 
ting the significance level for the one-tailed Spearman-Rank 
test at 99%, we find that under model (i) only 0.1% of sim- 
ulated samples could produce such a result, whereas under 
model (ii) approximately 94% of samples can do so. Model 

(i) may be rejected at a high level of significance, but model 

(ii) is consistent with the data. The physical interpretation 
of this result is unclear. 



3.2 Jets 

We begin by treating the upper limits on jet prominence as 
measurements and considering only beaming to be respon- 
sible for the distribution of observed jet prominences. We 
take the spectral index of a jet to be 0.8; the analysis is not 
significantly affected if this number is altered. The results of 
the simulations are plotted in Fig. |^. There is a clear best- 
fit value, corresponding to /3j = 0.62, pj = 0.005, although 
again a wide range of parameters is formally acceptable. As- 
signing a confidence range in the standard way, we find that 
/3j = 0.62j|g'j'* at the 90 per cent confidence level. For a 
(^4^nified model with critical angle 45°, we would expect the 
median BLRG jet with this value of /3-, to be about 4 times 
more prominent than the median NLRG jet; this appears 
consistent with the distribution seen in Fig. ^, though it is 
somewhat higher than the ratio (~ 2) of the medians quoted 
in paper I, where the median jet prominence was calculated 
on the assumption that all the upper limits fell below the 
median. The data are better fitted by smaller values of /3j 
and/or larger critical angles. 

If we were able to replace our upper limits with true 
measurements, we would expect the distribution of jet 
prominences to broaden, so that the best-fit values of /3j 
would increase. Although we do not know how far below the 
upper limits the true measurements lie, we can assess the ef- 
fect that they have on the results of the fits by allowing them 
to vary while the fits are being performed. We therefore ran 
a set of fits in which, at the same time as the comparison 
population was generated, each of the 7 upper limits present 
in the data was divided by some factor randomly drawn from 
a uniform distribution between 1 and 50. The value 50 was 
chosen as being the ratio between the highest jet prominence 
upper limit and the lowest actual measurement, though we 
expect that the difference between the upper limits and the 
true measurements will usually be less than this. The effect 
is to increase the best-fit values of /3j , but only slightly; this 
shows that the effect of the upper limits on our results is 
not a serious one. 

Whereas there is evidence from superluminal motions 
that relativistic beaming should dominate the prominences 
of cores, we have little direct evidence that it dominates the 
prominences of the kpc-scale jets, as discussed in section |l| It 
is thus worth considering what effect random variation of in- 
trinsic jet prominence (due to 'weather', as discussed above) 
might have on the best-fit beaming speeds. In general this 
will reduce the best-fit /3j, because less beaming is required 
to produce the observed spread in prominences; clearly if 
the distribution of intrinsic prominence is made to match 
the observed prominence distribution then no beaming is 
necessary. A further constraint on prominence distributions 
is provided, however, by the distribution of jet sidednesses 
in the data; the sidedness of a source is defined, as in pa- 
per I, as the ratio of the jet and counterjet prominences. 
If we assume that the intrinsic prominences of the jet and 
counterjet are drawn from the same distribution, then we 
can determine the distribution of sidednesses in a simulated 
dataset. Because we only have a few counterjet detections in 
the real dataset most of the sidedness data we have come in 
the form of lower limits, while we have no constraints at all 
on the sidedness of the 7 objects where no jet was detected. 
We therefore incorporate the sidedness constraint by mea- 
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suring the fraction of the simulated sidednesses that match 
the measured ones, either by being within ±20% of them if 
they are actual measurements or by being greater than them 
if they are upper limits. 7 of the simulated sidednesses are 
allowed to take any value. We then multiply the K-S prob- 
abilities of a fit to the prominence data by this fraction to 
give a (notional) total probability of matching the data. It 
turns out that for Pj ^ 0.5 almost all the simulated datasets 
can match the real dataset's sidedness in this sense; that is, 
the weighting fraction is 1 for /3j J; 0.5 and falls off below it. 

Unsurprisingly, then, when we draw the logs of the 
simulated intrinsic prominences from a normal distribution, 
which is (as can be seen from Fig. ^) not badly matched to 
the data, allowing both the mean (pj) and standard devia- 
tion (a) of the distribution to vary, the best fits are achieved 
with f3j ~ 0.5, a ~ 0.3 (logs being taken to base 10). How- 
ever, acceptable fits (probabilities > 0.05) are found for 
cr ~ 0.4, Pj ^ 0.1. We cannot use the data on jet promi- 
nences and sidednesses alone to rule out a model in which 
intrinsic jet and counterjet prominences are drawn at ran- 
dom from a sufficiently broad distribution with no or negli- 
gible relativistic beaming effects. If, however, we make the 
(perhaps more physically reasonable) stipulation that intrin- 
sic jet and counterjet prominences are the same in particular 
sources, then the simulations favour the relativistic beaming 
model more strongly, and the best-fit I3j is higher. 

We note finally that it is likely that jets do not emit ex- 
actly isotropically in their rest frames; this would affect the 
beaming relation. Where good polarization measurements 
are available the magnetic field in FRII jets appears to be 
directed longitudinally, though it is likely, given that the 
polarizations of jets are considerably less than the limiting 
value of ~ 70 per cent, that we are actually seeing a pro- 
jected, tangled field. In the limiting case of a uniform longi- 
tudinal field the correct prominence relation for jets is given 
by 

P; =P.[7.(l-/3.cose)]-(3+^°)(sine)i+" (5) 

(e.g. Cawthorne 1991), producing a different expected distri- 
bution of prominences for a given value of /3j. The increase 
in the exponent of the Doppler factor and the additional 
term in sin 9 to some extent cancel each other's effects, and 
so the cumulative probability distributions are not actually 
very different over the whole possible range of 9 for the mod- 
erate values of (3 that seem to be appropriate here. When we 
performed our fits using the prominence relation of equation 
BI we obtained very similar best-fit values of /3j but slightly 
higher values of pj . If the magnetic field is not uniform and 
longitudinal, as expected, then our analysis will be even less 
strongly affected. 



3.3 Relations 

In Fig. ^ we plot the core and jet prominences for the 31 
high-excitation objects in the sample. A correlation is ap- 
parent, significant at the 95 per cent level on a modified 
Spearman Rank test (as implemented in ASURv). As we ar- 
gued in paper I, if the prominences of cores are dominated 
by relativistic beaming, such a correlation is not expected 
unless jets too are beamed; any intrinsic relationship be- 
tween core and jet would be washed out by the beaming 
of the cores. But there is clearly substantial deviation from 



the predictions of the simplest relativistic beaming models, 
where intrinsic jet and core prominence are the same for all 
sources; in such models the points in Fig. |^ would be ex- 
pected to lie on a single curved line (cf. Bridle et al. 1994) 
which we plot, for comparison, in Fig. ^ Scatter about this 
line can arise in a number of ways. It is possible, for exam- 
ple, that there is normally jet bending on 100-pc scales, so 
that the beaming angle 9 is different for core and jet; it is 
also possible that /3c and l3j vary significantly from source 
to source. In the simplest models, though, the scatter about 
the solid line in Fig. |^ has implications for the distribution 
of intrinsic core and jet prominences, discussed above. 

We have already argued that the cores are necessarily 
beamed, but that the data seem to show a relationship be- 
tween the intrinsic core prominence and the linear size of 
the source — a relationship which is not obvious in the jet 
prominences (see Fig. 5 of paper I). We have seen that jet 
prominences are well modelled with a combination of intrin- 
sic scatter and beaming. We may now ask whether these two 
models, when combined, are consistent with the observed 
relationship between core and jet; and we can try to make 
quantitative our conclusion that models without beaming in 
the jets are ruled out by this relationship. 

We begin by asking whether the combined models are 
consistent with the observed correlation found on the Spear- 
man Rank test. Spearman's r for the 31 real data points, as 
calculated taking limits into account with ASURV, is 0.340. 
We simulate core and jet prominences based on equations ^ 
and 1^ respectively, drawing the logs of the intrinsic promi- 
nences of simulated jets from a normal distribution as in 
section and drawing the lengths of radio sources from 
the distribution discussed in section for cores, as be- 
fore, we choose f3c = 0.98. Unsurprisingly, models in which 
the spread in intrinsic prominence dominates the distribu- 
tion of jet prominences (so that jet and core are essentially 
unrelated) cannot adequately reproduce the correlation; if 
Pj = 0.1, a — 0.4, the probability of obtaining a correlation 
at least as strong as this is ~ 5 per cent, the level expected 
to occur by chance, whereas if Pj — 0.5, a — 0.3 it is a more 
acceptable 75 per cent. We can rule out at the 90 per cent 
confidence level all the best-fit models for jet prominence 
derived in section 3.2 that have Pj ^ 0.2, reinforcing our 



conclusion that beaming is necessary in the jets. 

As a final test, we can ask whether the distribution 
of core and jet prominences in Fig. ^ is consistent with 
the model we have devised. Peacock (1983) provides a two- 
dimensional version of the Kolmogorov-Smirnov test, subse- 
quently refined by Fasano & Franchesini (1987); we use the 
implementation of their test given by Press et al. (1992), 
though we have verified that the original version of Peacock 
(1983) gives similar results for our data. Fixing the parame- 
ters for the core at Pc = 0.98, pc ~ 0.05, we allow the jet Pj 
and intrinsic prominence to vary; we find, as we expect, that 
the results are similar to those plotted in Fig. ^ The data 
are (just) acceptably fitted by a model with no intrinsic scat- 
ter, corresponding to the solid line in Fig. ^. However, they 
are better fitted (K-S probabilities of 50-60 rather than 10 
per cent) by models in which there is some dispersion in the 
intrinsic prominences. If we allow the intrinsic core promi- 
nences to depend on length, as above, and the intrinsic jet 
prominences to be drawn from a normal distribution, then as 
before the best results are obtained with Pj ~ 0.5, a ~ 0.3. 
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Models in which the jet and core are uncorrelated {6c and 
9j are both randomly chosen rather than being the same) 
give acceptable but poorer fits to the data, while models in 
which the jet and core are anticorrelated (9c = (7i"/2) — 9j) 
are never acceptable fits, as we would expect. 



4 DISCUSSION AND CONCLUSIONS 

The results of the previous section allow us to conclude with 
little ambiguity that relativistic beaming really is signifi- 
cant in determining the prominences of the kpc-scale jets of 
these objects, providing that it also dominates the promi- 
nences of the cores. Speeds between 0.5 and 0.7c seem to 
provide the best fits to the data. Speeds much higher than 
this do not reproduce the distribution of jet prominences 
well, while lower speeds cannot well reproduce the observed 
relationships between core and jet prominence and the jet 
sidedness distribution, and are in any event at odds with 
observations of the Laing-Garrington effect. These conclu- 
sions are similar to those in earlier work (Bridle et al. 1994; 
Ifough 1994; Wardle & Aaron 1997) but they are based on 
a sample that should be free of (unknown) orientation bias. 
Subject to the caveats already discussed, the data and the 
relativistic beaming speeds derived from them are broadly 
consistent with the simple unified model of paper I, in which 
BLRG are aligned at angles less than ~ 45° to the line of 
sight. 

Are these beaming speeds a reflection of any real phys- 
ical velocity in the beams? As pointed out by Bridle et al. 
(1994), it is entirely possible that the material in the beam 
does not move with a single speed. If this is true, the speed 
derived from a single object is, as they say, a weighted av- 
erage over the velocity field in the beam; the speed derived 
by application of equation ^ to a whole sample of beamed 
objects is still more complicated in its relationship to any 
true physical speeds. Related to this is the whole question 
of where the jet emission comes from, and why, in fact, we 
see emission from the beam at all. However, if (as has been 
suggested; e.g. Bridle 1996) the emission in these objects 
comes from a shear layer at the edge of the beam, while the 
emission from a higher-velocity central spine is suppressed to 
negligible levels by beaming, then the characteristic speeds 
we derive may well be a reasonable estimator of the physi- 
cal speeds in the shear layer. Such a picture may also help 
us to understand why the speeds derived from cores and 
jets are so different; perhaps we are simply seeing different 
velocity regimes because different regions of the beam are 
emitting. Alternatively, it may be that the difference is evi- 
dence of a genuine deceleration of the beam between the two 
scales. Such deceleration is invoked in models of FRI radio 
sources, where similar speeds around 0.5c are derived, from 
sidedness arguments, at the bases of the jets; [e.g. Hard- 
castle et al. (1996), Laing (1996), Hardcastle et al. (1997)]. 
As discussed by Bowman, Leahy & Komissarov (1996) it is 
possible to decelerate FRII beams by entrainment without 
excessive dissipation. 

Another, related question is raised by the derived intrin- 
sic prominences. The best-fit intrinsic prominences of cores 
in our sample are ^ 0.02, with the value for our represen- 
tative Pc ~ 0.98 (7c ~ 5) being about 0.05. These numbers 
are similar to those that have been derived in similar ways 



from other samples [e.g. Morganti et al. (1995) and refer- 
ences therein]. But the best-fit intrinsic prominences of jets 
are lower; pj ~ 0.005 is a typical value derived from the 
fits for /3j ~ 0.6. In other words, even if the sources were 
not beamed, the sub-kpc-scale radio emission from the cores 
would be brighter than that from the much larger 100-kpc- 
scale jets. The beam is evidently much less efficient during 
the first few hundred pc of its journey to the hot spots than 
it is at any future time. This may be a sign either of strong 
interaction with the external medium or of shocks in the 
jet on these small scales, which would be consistent with 
a model in which the beam is decelerating strongly at this 
stage. 

It is worth noting that the efficiency of the jet in these 
FRIIs is still very high; because the total radio luminosity 
of an extended source is much lower than the total power 
supplied by the beam (cf. Rawlings & Saunders 1991) and 
the unbeamed luminosity of the core and jet is less, by ap- 
proximately Pj +Pc, than the total luminosity of the source, 
it follows that the efficiency of the jet is [1 — {pj + Pc)], 
or probably > 99 per cent. Almost all the energy supplied 
by the beam is used to excite the lobe electron population 
and to do work on the external medium, and very little of 
it is radiated away in jets. 
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Table 1. Fluxes and prominences for the sample sources 
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Column 2 gives the observing frequency of the radio data used in the analysis. Prominences are corrected to a nominal frequency of 8.4 
GHz as described in paper 1. Column 3 lists the emission line type of the source. 'E' indicates a LERG, 'N' a NLRG and 'B' a BLRG; 
'?' indicates an unclassified source. Paper I provides the references for these classifications. Column 4 gives the total flux density at the 
observing frequency. Errors on this quantity (column 5) axe quoted only where the flux density was taken from the literature rather than 
VLA maps, or where there was some uncertainty in the definition of the source region on the VLA maps, as discussed in paper I; the 
standard VLA calibration errors apply to all other measured quantities. 
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Figure 1. Core prominences of the objects in the sample classed by emission-line type. Dashed shading indicates an upper limit. 
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Figure 2. Jet prominences of the objects in the sample classed by emission-line type. Dashed shading indicates an upper limit. 
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Figure 3. The relationship between core and jet prominences for the high-excitation objects in the sample. Filled symbols are BLRG, 
others are NLRG; arrows indicate upper limits. The solid line shows a th eore tical curve assuming only relativistic beaming affects the 



prominences, using f3j = 0.62, pj = 0.005, /3c = 0.98, pc = 0.05; see section 3.3 
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Figure 4. Kolmogorov-Smirnov probability contours for model core prominence data fitted to the prominences of BLRG and NLRG. 
Beaming speed /9c is plotted on the x-axis. Intrinsic prominence pc is plotted on the y-axis. Contours at K-S probabilities of 
0.05, 0.1, 0.15 ... , 0.40. Above, /3c ranges from 0.1 to 0.98; below, as above but showing the region with 0.92 < f3c < 0.998. 
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Figure 5. Kolmogorov-Smirnov probability contours for model core prominence data fitted to the prominences of BLRG and NLRG, 
including a dependence of intrinsic prominence on source length as discussed in the text. Beaming speed 3c is plotted on the a;-axis, and 
varies between 0.92 and 0.998; the constant relating length to intrinsic prominence, fc<;,i (kpc ~^), is plotted on the j/-axis. Contours at 
K-S probabilities of 0.05, 0.1, 0.15 ... , 0.75. 




Figure 6. Kolmogorov-Smirnov probability contours for model jet prominence data fitted to the prominences of BLRG and NLRG. 
Beaming speed /3j is plotted on the x-axis, and varies between 0.1 and 0.98; intrinsic prominence pj is plotted on the y-axis. Contours 
at K-S probabilities of 0.05, 0.1, 0.15 ... , 0.70. 
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